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Abstract. We show that sudden variations in the composition and structure of an hybrid star can be triggered by its rapid 
spin-down, induced by r-mode instabilities. The discontinuity of this process is due to the surface tension between hadronic and 
quark matter and in particular to the overpressure needed to nucleate new structures of quark matter in the mixed phase. The 
consequent mini-collapses in the star can produce highly energetic gravitational wave bursts. The possible connection between 
the predictions of this model and the burst signal found by EXPLORER and NAUTILUS detectors during the year 2001 is also 
investigated. 
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1. Introduction 

The analysis of the data collected by the Gravitational 
Wave antennas EXPLORER and NAUTILUS during the year 
2001 shows the existence of coincidences between the sig- 
nal detected by these two resonant bars ( |Astone et al. (20 02) 
|Coccia et al. (2004)) . In particular, an excess of coincidences 
respect to the background is concentrated around sidereal 
hour four, which corresponds to the orientation for which 
the sensitivity of the bars is maximal for a signal coming 
from the direction of the galactic center. Although the statis- 
tical significance of this signal is debated (see |Finn (2 003) 
|Astone et al. (2003b)||Astone et al. (2003a)) , it is interesting to 
investigate the possible origin of the inferred signal using ex- 
isting models of GW emission. In particular we will concen- 
trate on the ones associated with instability modes of a rotat- 
ing compact star. It is important to remark that in the analy- 
sis performed on the data from the two gravitational antennas 
only impulsive signals have been searched for, having a typ- 
ical duration of ~ 10 ms. In the literature an extensive anal- 
ysis has been performed concerning steady sources and peri- 
odic GW emission (|Andersson et al. (2002)| |Wagoner (2002)| 
|Reisenegger & Bon acic (2003)). In this paper we will explore 
the possibility that a transient GW burst is triggered by the 
steady emission of GWs generated by r-mode instabilities. In 
particular, we will show that the almost continuous dragging 
of angular momentum from the star can induce a sudden vari- 
ation in its structure and composition, generating a few bursts 
of GWs. 



The first ingredient of our model are the so called 
r-mode instabilities, which are very efficient in drag- 
ging angular momentum from a rapidly rotating compact 
star ( jAndersson (1998)| [Friedman & Morsink (1998)) . A rapid 
spin-down deeply affects the structure and composition of 
the star. This problem has been discussed in detail in a 
few papers in the past years (see |Glendenning et al. (1997)| 
|Chubarian et al. (2000)} , where it has been shown that, using 
rather standard values for the model parameters, the composi- 
tion of an hybrid hadronic-quark star can be dramatically mod- 
ified when the rotational frequency changes in a range centered 
around a few hundreds Hertz. In particular, during the spin- 
down era the central density increases and a larger amount of 
matter is converted into a mixed hadron-quark phase (MP) or 
into pure quark matter. The new idea we discuss in this pa- 
per is the effect of a surface tension at the interface between 
hadronic and quark matter. If the surface tension is not van- 
ishing, the formation of new structures in the MP can be de- 
layed. A similar idea has been discussed recently in connec- 
tion with the transition from a metastable hadronic star into 
a stable hybrid or quark star, showing that a huge amount of 
energy can be released in the process of conversion of bulk 
hadronic matter into quark matter (jBerezhia ni et al. (2 003) 
Prago et al. (2004)| |Bombaci et al. (2004)) . The mechanism 
discussed in the present paper is also based on metastability, 



1 A similar idea has recently been discussed in 
Ref. jHarko et al. (2 004) l. There, the possibility that the spin 
down can trigger a phase transition is analyzed taking into account the 
effect of the surface tension. In Ref. I Hark o et al. (200 4) l the angular 
momentum of the star is dragged by the magnetic dipole radiation 
which can be an efficient mechanism of spin down only in presence 
of a strong magnetic field. 
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but it involves phase transitions in the MP, which are respon- 
sible for a relatively large modification of the stellar structure 
without a large energy release. We suggest the possibility that 
the formation of the MP can take place in several steps, in each 
of which the radius and the oblateness of the star change only 
by a relatively small amount. A crucial observation is that the 
conversion process from hadronic matter into MP will start in 
one point inside a spheroidal layer of metastable matter (to be 
discussed later) and will propagate inside that same layer at fi- 
nite velocity. In this way the burning process generates a large 
non-radial oscillation of the star. During this period a strong 
GW emission takes place, and we will discuss its possible con- 
nection with the analysis of the gravitational bars data by the 
ROG Collaboration 2 . 



2. R-mode induced spin down 

In the literature, two scenarios have been discussed concerning 
a steady GW emission due to r-mode instabilities. The first sce- 
nario is based on the emission of GWs from a hot and rapidly 
spinning compact stellar object, which has not lost its angu- 
lar momentum in the very first part of its existence after the 
supernova explosion. This is possible if the bulk viscosity is 
large for temperatures of order ~ 10 9 K or higher, so that r- 
modes are damped till the temperature drops below ~ 10 9 K. 
Possible candidate s are quark stars (jA ndersso n et al. (20 02) I, 
hyperonic stars ( Lindblom & Owen (2002) ) or hybrid stars 
Prago et al. (2005)| l. The second scenario involves older stars 
which are reaccelerating due to mass accretion from a com- 
panion. In this way a sort of "cycle" can develop (see Fig.Q, 
in which the star periodically goes through the following 
steps: 1) mass accretion with increase of angular velocity; 
2) instability due to r-mode excitations with reheating due 
to bulk viscosity; 3) loss of angular momentum with emis- 
sion of GWs. The reheating due to bulk viscosity during 
phase 2) is so efficient that, if the instability region is reached 
from the low-temperature side, the star rapidly reheats and 
reaches the high-temperature side of the instability region 
( Reisenegger & Bonacic (2003)| l 3 . 

Let us introduce the equations regulating the process 
of GW emission in the two scenarios discussed above. 
Here we follow the formalism of Wagoner (2002) and 
Kinney & Mendell (2003) where the scheme proposed by 
OwerTeUd"*(T998')|has been refined 4 . The evolution equations 



2 Another possible relation between phase transitions and GW 
emissions is discussed in Ref. (Mi niutti et al. (2003)) . In that paper 
phase transitions, associated with density discontinuities, can excite 
the so called g-modes and large non-radial oscillations can develop. 

3 Also neutron stars can recycle, as suggested in 
Refs. (Levin (1999)1 jAndersson et al. (2000)||Heyl (2002)) . 

4 We assume that the r-mode canonical angular momentum con- 
tributes zero physical angular momentum to the star. This is 
a sensible approximation as results both from explicit calcula- 
tions (Levin & Ushomirsky (20017} and also from the observa- 
tion that the r-modes saturate quickly and after saturation no 
physical angular momentum can go into them fWagoner"(2 002) 
IKinney & Mendell (2003)) . 
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Fig. 1. r-modes instability window (upper panel). Coming 
from high temperatures (A) the star loses its angular mo- 
mentum until it exits the instability region (B) (first sce- 
nario). By mass accretion the star can be re-accelerated 
until it reaches the instability window from the low tem- 
perature side (C). The excitation of r-modes instability 
produces a fast reheating of the star due to bulk viscos- 
ity dissipation (D) (second scenario). In the lower panel a 
magnification of the regions indicated by (A) and (D) is 
displayed. Here we used a mass accretion rate of 1O~ 8 M 
per year. 



read: 




Here a is the dimensionless amplitude of the r-mode, l/t v = 
l/t s + l/tb, t g , t s and % are time scales associated with GW 
emission, to shear and to bulk viscosity damping, respec- 
tively. / and J are dimensionless values of the moment of in- 
ertia and of the angular momentum (for all details see e.g. 
Ref. ( Andersson et al. (2002)i). Eq. (1) describes the damp- 
ing of r-modes due to viscosity, Eq. (2) describes the varia- 
tion of the angular momentum of the star in presence of the 
torque due to gravitational wave emission and to mass accre- 
tion and, finally, Eq. (3) describes the thermal evolution given 
by the contributions of the reheating due to mass accretion, 
shear and bulk viscous dissipation of the r-modes and of the 
cooling due to neutrino emission. Obviously, in the first see- 
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nario depicted above, mass accretion is not present. To com- 
pute the time scale associated with bulk viscosity for hybrid 
stars we use the results of Ref . ( Dra go et al. (2005)) in which 
the viscosity of MP has been computed. In that paper it has 
been shown that the viscosity of MP is of the same order of 
magnitude as the viscosity of pure quark matter if supercon- 
ducting gaps are not present or it is reduced by a factor ~ 10 
if a color superconducting 2SC gap is taken into account (see 
Madsen (2000) I. Concerning the value of shear viscosity, we 
have taken into account not only the contribution associated 
with pure quark matter (Anders son et al. (2002)) , but also the 
contribution associated with the viscous boundary layer which 
is present in a star having a crust made of nucleonic matter 
(Bildsten & Ushomirsky (1999)1. In Fig.[2]we plot the dimen- 
sionless amplitude h and the frequency / of the GWs emitted 
by a a star in the first scenario. In Fig. [3] the same quantities 
are displayed for the second scenario. In this figure, a mass 
accretion rate of order 10~ 8 M o per year has been assumed. If 
the value of the accretion rate is reduced, the star re-enters the 
instability window at smaller temperatures and larger values 
of the angular velocity. In this way, the GW emission in the 
second scenario would be rather similar to the one obtained in 
the first scenario. Finally, let us remark that the same compact 
star can enter first the instability window associated with the 
first scenario, when the star is relatively young, and, after some 
time, the star becomes again unstable due to mass accretion as 
described in the second scenario. These possibilities are dis- 
played in Fig.^ 

As it can be seen in Figs.[2j|3] in both scenarios the initial 
part of the signal decomposes into bursts lasting a few minutes 
and separated by periods of a few days (or months) of quies- 
cence. During this initial phase of the emission, the star fol- 
lows a trajectory in the temperature - angular velocity plane, 
oscillating around the instability line displayed in Fig. [2 (see 
also Ref. (Anders son et al. (2002)) ). After this phase, which 
can last months or years, the angular momentum is dragged 
almost continuously and the signal becomes steady for hun- 
dreds of years until the star finally exits the instability region. 
In principle, the first part of the signal having an amplitude h 
~ 10~ 22 for a source located at 1 kpc, could be detected by 
resonant bars. The main difficulty associated with the search 
of this type of periodic signal is due to the almost monocro- 
maticity of the GW emission which could be detected only 
by a dedicated search. For instance, in the case of the pul- 
sar J1939+2134, a dedicated search has been performed for a 
signal having a frequency twice the rotation frequency of the 
star (which would correspond to the signal emitted by a stel- 
lar object having a non-axisimmetric shape), with a null result 
(see |Abbott et al. (2004)) . On the other hand, the signal asso- 
ciated with r-modes would mainly be emitted at a frequency 
f g = m/3n. 

3. Toy-model for a rotating star 

The r-mode mechanism we have discussed so far produces a 
steady emission of periodic GWs and therefore it does not cor- 
respond to the signal discussed in Refs. ( |Astone et al. (20 02) 
|Coccia et al. (2004)) . The new idea we are now introduc- 
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Fig. 2. The frequency (upper panel) and amplitude (lower 
panel) of the GW signal emitted by r-mode instabilities 
as functions of time. Here the compact star is approach- 
ing the instability window from the high temperature side 
(first scenario). The computed amplitude corresponds to 
a distance of 1 kpc. 



ing is a possible connection between the rapid dragging of 
the angular momentum and discontinuous changes of the 
structure of the star. It is worth remarking that r-modes 
are a very efficient way of reducing the angular momentum 
of a compact star. Although it will not be discussed here, 
also a strong magnetic field as the one present in magne- 
tars can efficiently dragg the angular momentum of the star 
by dipole magnetic radiation (j Harko eUjl^OOJQ }. As shown 
in Refs. ( |Glendenning et al. (1997)| |Chubarian et al. (2000)) , 
when a hybrid star slows down its density increases till the 
first critical density is reached at the center of the star and 
MP starts being produced. A further reduction of the an- 
gular momentum allows the star to increase the fraction of 
MP occupying its core, until the central density reaches the 
second critical value. At that point pure quark matter starts 
being produced. During all the spin-down process the ra- 
dius of the star decreases by several kilometers. This gradual 
modification of the structure of the star can take place dis- 
continuously in several (small) steps if the effect of a non- 
vanishing surface tension is taken into account. The inner 
regions of the star in which MP will later form becomes 
firstly metastable. The formation of germs of stable MP pro- 
ceeds through quantum tunneling. The probability of this 
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Fig. 3. Same as in fig. 1. Here the compact star is ap- 
proaching the instability window from the low tempera- 
ture side after the " recycling" process (second scenario). 



transition strongly depends on the value of overpressure. As 
shown in Refs. (|Iida & Sato (1998)| |Berezhiani et al. (2003)| 
|Drago et al. (2004)||Bombaci et al. (2004)1 the nucleation time 
needed to form new structures of quark matter can be very 
long if the overpressure is not large. When the overpressure 
reaches the critical value, in one randomly chosen site inside 
the metastable layer a new drop of quark matter forms. As it 
will be shown, the process of conversion of hadronic into quark 
matter propagates with finite velocity v c inside the star and a 
sudden modification of the composition and of the structure of 
the star occurs during a timescale ~ R/v c . During this period 
non-radial modes develops and a few bursts of GWs can be 
emitted until a new equilibrium configuration is reached. 

In order to give a qualitative estimate of the magni- 
tude of the oscillations and therefore of the amplitude of 
the GWs emitted we resort to the toy model proposed in 
Ref. Peiselberg & Hjorth-Jensen (1998)| l. We model the hy- 
brid star as a spheroid containing a MP core with uniform 
density p2 and a crust of nuclear matter with uniform density 
pi. The pressure at which the phase transition from the first 
to the second component takes place is a parameter Pq. The 
Newtonian hydrostatic equation reads: 



1 dP m(d) 2 2 
-— = -G — — + -O a, 

p da a L 5 



where 



m(a) — An 



Jp(a')a' 2 da . 
o 



(5) 



Here the "effective radius" a appears, which is related to a po- 
sition inside the star by the equation: 



r = a(l - eP 2 (cos6)), 



(6) 



where e = 



5 n 2 



(4) 



-j w t and P2(cos6) is the 2nd Legendre poly- 
nomial. Eq. (4) is then solved and in this way the pressure 
inside the star is analytically determined for a given angu- 
lar velocity. In Fig. 0] the structure of a rotating star is dis- 
played. Comparing with results based on the actual solution 
of the equilibrium equation for a rotating compact object 
(see e.g. jGlendenning et al. (1 997)1 |Chubarian et al. (2000)> ), 
it is possible to notice that the toy model is at least qualita- 
tively in agreement with the exact solution. Unfortunately, the 
toy model badly underestimates the quantitative effect of the 
spin-down on the radius and composition of the star. For in- 
stance, in the toy model we get reductions of the radius of 
order 0.5 km, while the calculations based on realistic equa- 
tions of state give reductions typically of order a few km 
UGlendenning et al. (1 997)1 |Chubarian et al. (2000)> . The same 
problem arises for the estimate of the fraction of the volume of 
the star occupied by MP: the toy model gives an increase of this 
quantity of roughly 20%, while more realistic calculations give 
a significantly larger value. It is clear that we are at the moment 
only aiming at a qualitative study of the proposed mechanism 
and that refined calculations are needed in order to get a more 
precise estimate. 

The question we want to numerically investigate within the 
toy model is the following: given a specific equation of state 
(in our case fixing the values of pi, pi and Pq) and a value 
for the surface tension <x, which variation of the angular veloc- 
ity AQ./Q is large enough to trigger the formation of a critical 
drop of quark matter, in a time scale of order days or years? The 
crucial ingredient in this calculation is the relation between the 
overpressure AP/P and AQ./Q.. The overpressure is determined 
computing, for a same element of fluid, the difference between 
the value of the pressure after and before the slow down (la- 
grangian perturbation). In particular, we are interested in the 
value of the overpressure in the region immediately surround- 
ing the core of already formed MP. In Fig.[5]we show the over- 
pressure, in the above defined region, for a value AQ/Q = 0.05 
and for various values of O. As it appears, AP/P ~ a AQ/f2, 
with a ~ 0.3 4- 1.4, and the larger values correspond to the 
faster rotating stars. It is also interesting to notice that, as the 
star slows down the layer in which new MP is formed moves 
to outer regions and the thickness of the layer shrinks. In Fig. [6] 
the overpressure inside the star is displayed for a fixed value of 
£2 and AQ/£L Notice that the value of the overpressure is larger 
in the layer around the MP core and therefore in that region the 
conversion will take place with a larger probability. 

We need now to compute the nucleation time for the 
obtained value of the overpressure and a given value of 
cr. This can be done following the formalism developed in 
Refs. ( |Iida & Sato (1998)| |Berezhiani et al. (2003)| > and based 
on quantum tunneling. Within the toy model, nucleation 
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Fig. 4. The total radius of the star and the radius of the 
inner core of MP are shown as functions of the angular 
velocity in the toy model. 
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Fig. 5. The value of the overpressure as a function of the 
effective radius, for four values of the angular velocity. 
Here AQ/Q = 0.05. 



4. Emission of GW bursts 

4.1. Time structure of the bursts 

We can now describe more precisely the GW emission in the 
mechanism here introduced. As it has been shown in Figs.fJJ- 
[5] it is possible to reduce the angular velocity of the star by 
some 10-20% in ~ 10 years via emission of periodic GWs in- 
duced by r-mode instabilities. This reduction corresponds to an 
increase of the inner pressure by roughly the same percentual 
amount. As we have seen, when an overpressure of order a 
few percent is reached, the star will reassess forming a new 
region of MP. Therefore, we can expect a few periods of GW 
bursts activity in 10 years. It is interesting to remark that in 
this model the mechanism of emission of GWs is similar to 
the one invoked in connection with the so-called soft-gamma 
repeaters (see de Freitas Pacheco (1998) ). In particular, in both 
models the burst of GWs is due to a starquake activity. The 
main difference between the two models is that, in the soft- 
gamma repeaters case, the origin of the sudden change of the 
structure of the star is related to the maximum shear stress that 
the crust can bear before cracking. In the model proposed in 
the present paper, the role played by the rigidity of the crust 
is now played by the surface tension at the interface between 
hadronic and quark matter. As we will show, it is possible in our 
model to obtain more violent oscillations of the star and there- 
fore more energetic GW bursts. This is not surprising since the 
physics underlying the soft-gamma repeaters model is based 
on a metastability related to the typical energy scale of atomic- 
nuclear physics, while in our model the energy scale is related 
to hadronic physics. 

In the model here described the emission of bursts of GWs 
is due to the collapse of the region involved in the formation of 
new MP. It is not easy to estimate the exact size of the interested 
region. Once a drop of quark matter develops in the zone where 
MP is being produced, the neighboring areas are probably trig- 
gered so that a relatively large fraction of the metastable layer 
will collapse. We can imagine that, after the mini-collapse, a 
compression wave develops and propagates from the site in 
which the new quark matter drop formed to the neighboring re- 



time of order days can be obtained, for the values of 
the overpressure as the one displayed in Fig. [5] and us- 
ing values of cr ~ a few MeV/fm 2 . These values of <x 
are of the same order of the one estimated in the MIT 
bag model (Berg er & Jaffe (1987)) . They are also not far 
from the values investigated in Refs. |Berezhia ni et al. (2003)| 
|Drago et al. (2004)| |Bombaci et al. (2004)) . We would like to 
stress again that refined calculations are needed in order to get 
a quantitative estimate of the nucleation time. The most arti- 
ficial feature of the toy model we are using is probably re- 
lated to the absence of modifications of the already formed MP 
during the spin-down process. In a realistic model, not only 
the volume occupied by MP increases during the spin-down, 
but also new structures form, with an evolution of the existing 
structures from drops, to rods, to slabs as the star slows down 
( |Glendenning (1997)| . 



Y [km] 




Fig. 6. The value of the overpressure inside the star, for 
Q = 5000 s" 1 and AQ/Q = 0.05. Lighter areas corre- 
spond to higher values of the overpressure. 
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gions as shown in the picture of Fig. The compression wave 
corresponds to a temporary overpressure that adds on the al- 
ready existing pressure inside the metastable matter. To inves- 
tigate if the propagation of the compression wave is able to trig- 
ger the formation of new structures of quark matter, we com- 
pare two timescales. The first timescale f trans corresponds to the 
duration of the overpressure Sp, due to the compression wave, 
on a specific point inside the metastable layer. It can be com- 
puted from the velocity of the compression wave v c (which is 
the velocity of sound and it should be a fraction of the velocity 
of light) and the size of the compression wave itself A, which 
should be of the same order as the radius of the newly formed 
quark matter drop whose formation originated the wave, i.e. 
A ~ 10 fm. In this way we get f trans ~ 10~ 20 s, where only an es- 
timate of the order of magnitude is needed. This timescale must 
be compared with the nucleation time f nuc i that corresponds to 
the value of overpressure obtained adding the already existing 
overpressure in the metastable layer Ap to the overpressure due 
to the compression wave 5p. Assuming dp ~ Ap and using the 
standard formalism we obtain t mc \ < ftrans- We conclude there- 
fore that the conversion process can indeed propagate inside 
the metastable layer with a velocity marginally lower than v c . 

Within our model, it is not possible to give a precise esti- 
mate of the size of the region collapsing in a timescale of order 
1 s or smaller. The extension of this region can vary from the 
whole metastable layer to a relatively small portion of the star. 
For instance, it is possible that the conversion will propagate 
only in the area in which new drops of QM are not too far 
from each other. The randomness of the size of the collapsing 
region is typical of quake phenomena. Once a fraction of the 




Fig. 7. Picture of the mini-collapse process. The nucle- 
ation site of the new drop of MP is indicated by the small 
gray circle, v c indicates the velocity of the compression 
wave and 5v is the velocity of the infalling upper layers. 



metastable layer has collapsed, the other parts will presumably 
follow the same fate in a timescale much shorter than the time 
needed to reach the critical value of the overpressure. In other 
terms, we expect to have a few bursts taking place in a rela- 
tively short period, while a much longer delay (of order years) 
separates the phases of quake activity. This feature is similar 
to the temporal distribution observed in soft-gamma repeaters 
l |Cheng et al. ( 1995) I and interpreted there as due to starquake 
activity jdTFreitas Pacheco (1998)| l. Interestingly, the tempo- 
ral clustering of the events seems also to be a characteristic 
feature of the GW bursts emission jCoccia et al. (2004) I. A de- 
tailed analysis on the temporal structure of the GW bursters has 
been made in Ref. (Du bafh et al. (2003)| l where also the detec- 
tion strategies of these type of signals are discussed. 

4.2. Energy of the GW bursts 

We will now try to estimate the energy of the GW bursts. It is 
not easy to provide a realistic approximation to this quantity 
since it would require a detailed analysis of the dynamics of 
the micro-collapse. In this Section we will first present a rough 
estimate based essentially on dimensional analysis and, after, a 
more quantitative discussion which will provide a lower limit 
to the energy released. 

4.2.1 . Order of magnitude estimate 

The order of magnitude of the energy released can be estimated 
from the equation: 

E cw = M(AR/R) 2 , (7) 

where M is the mass in quadrupole motion and AR is the am- 
plitude of the oscillation. In our model, the entire region of 
the star above the layer where the conversion occurs would 
participate to the quadrupolar motion, since the external re- 
gions have to readjust immediately after the new layer of MP is 
formed. Therefore, we can assume that the mass in quadrupo- 
lar oscillation is of the same order of the total mass of the 
star. Concerning AR, it is of the order of the shrinking of the 
radius of the star due to formation of a new layer of MP. 
Using the toy model results, we can estimate that for each 
mini-collapse the variation of the radius of the star is of or- 
der 20-30 m, which corresponds to AR/R ~ 2 4- 3 x 10~ 3 . 
We can now put an upper limit to the energy released in GWs. 
Assuming that the oscillation induced by the phase transition is 
completely quadrupolar, the energy released would be of order 
Ecw ~ (0.5 t 1) x 1O~ 5 M . In the next Section we will study 
how the oscillation decomposes into the normal modes of the 
star and we will find that only a a few percent of the total en- 
ergy is associated with quadrupolar oscillations. On the other 
hand, a more realistic model should give larger variations of 
the radius for a same value of AO/Q. We can expect that in a 
more realistic calculation the upper limit to the energy released 
in GWs can be up to an order of magnitude larger, approaching 
E cw ~ 10 '/V/ : . . 
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4.2.2. Decomposition of the perturbation into 
nonradial modes 

In this Section we try to give a more quantitative estimate of 
the amplitude of the GW, by explicitly performing the decom- 
position of the perturbation into the normal nonradial modes of 
the star. In the following, we will first provide a "model" of the 
perturbation of the star's fluid induced by the conversion pro- 
cess. Then, after computing the spectrum of non-radial modes 
of the star, we decompose the perturbation on the basis of the 
normal oscillation modes, 77,,/. In particular, we can compute the 
f-mode component 77/2 of the perturbation and the correspond- 
ing GWs emission. This is important because the excitation of 
f-mode is the most efficient way for producing GWs. 

In order to model the perturbation of the pressure of the star 
during the micro-collapse, we resort again to the toy model 
and we assume that the perturbation corresponds to the over- 
pressure already calculated in Sec|5] since immediately after 
the conversion the overpressure reduces to zero. In Sec|3]we 
computed the overpressure AP/P as a function of the effective 
radius a, for fixed values of D. and AD./Q.. Inverting the relation 
given by Eq. we can write the overpressure as a function 
of r and 6. For small values of the deformation parameter e the 
perturbation can be approximated by: 



AP AP d(AP/p) 

p p ar 



(8) 



This corresponds to a decomposition of the perturbation 
into a radial term, which will not contribute to the emission 
of GWs, and a non-radial quadrupolar term, which will be fur- 
ther decomposed using the normal non-radial modes basis. In 
the following, while we compute the nonradial modes within 
a realistic model for the star, concerning the form of the per- 
turbation AP/p we assume that it is given by the expression 
obtained within the toy model. This crude approximation can 
be overcome by computing the structure of a rotating star us- 
ing a realistic EOS |Drago et a\.\ . A generic perturbation can 
be expanded as: 



The power emitted in GWs by the excitation of the f-mode 
can be estimated using the formula ( Ip ser & Lindblom (1991)| l: 



dE\ 
dt ) n 



-co 6 N 2 \AD 20 \ 2 



where the general expression for N/ reads 
G (/+!)(/ + 2) 



N[=4n 



c 2M 1(1- 1)[(2/+ 1)H] 2 ' 



(10) 



(11) 



a> is the frequency of the mode and ADi,„ is the mass multipole 
of the perturbed fluid. In our analysis we take into account only 
the quadrupolar contribution (I = 2). The mass multipoles read: 



AD 



Im 



-/ 



' Y* lm d x , 



(12) 



where Ap is the perturbation of the mass density. The function 
Ap can be expressed in terms of AP, using the relation 



p AP p AP 



(13) 



where F is the adiabatic index and c s is the sound velocity. 
Finally, for AP/p we use the / = 2, /-mode component ob- 
tained from eq.©. Therefore the density multipole reads: 



(Ap)/2 = ~iC fl T] fl Y ln 



(14) 



Substituting this expression in eqs. J1UI12I > we get an estimate 
of the power emitted in GWs due to the excitation of the f- 
mode: 



dE\ 

— -2.4x10 

dt ) f2 



_ 7 / oj/2tt\ 
\ 1kHz j 



M /s. 



(15) 



In our estimate of the power of GWs emission we have con- 
sidered up to now only the contribution of the mass quadrupole, 
but it is necessary to take into account also the contribution as- 
sociated with the current multipole, given by the expression: 



SJlr, 



(16) 



AP 



(r, Q,<f>)= > c„i m ri n i(r)Y Lm (6, </>) , 

P 1 



(9) 



where T]„i(r) are the eigenfunctions of the non-radial oscillation 
equation (see Ref. ( |Unno et al. (198 9))). The index n can indi- 
cate the fundamental / mode or the p and g modes. For simplic- 
ity, we have computed the spectrum of non-radial oscillations 
in the case of a non-rotating, Newtonian Hybrid Star, within the 
Cowling approximation and assuming that the Brunt- Vaisala 
frequency is vanishing. In this simplified case, g-modes have 
zero frequency, and the equation describing non-radial oscilla- 
tions reduces to a Sturm-Liouville form. 

In Fig.[S]we display the quadrupolar component of the per- 
turbation, defined in eq. (|8jl, together with the I = 2 eigenfunc- 
tions of the f-mode and of the first p-mode. It is now possible 
to compute the f-mode component of the perturbation, which 
turns out to be c/2 - 0.05. 



AP/p 



-0. 



-0.04 




Rfkml 



Fig. 8. The nonradial component of the perturbation AP/p 
(thick line) and the f-mode (continuous line) and first p-mode 
eigenfunctions (for 1=2), as functions of the radius of the star. 
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where Yf* are the magnetic spherical harmonics. In Eq. d!6t 
two terms are present: one, v 6p, is associated with the velocity 
of the fraction of the fluid which undergoes the density fluctu- 
ation . This term is suppressed respect to the one already esti- 
mated, due to the presence of a factor v/c which is obviously 
always smaller than one. The second term, p 6v{r), corresponds 
to the variation of velocity of all fluid elements inside the star. 
In our model, Sv(r) is associated with the rearrangement of the 
external parts of the star, immediately after the formation of 
a new layer of MP. It is very difficult to estimate Sv(r), and in 
particular its dependence on the position inside the star. We can 
assume that all regions above the conversion layer will move, 
as pictorially represented in Fig. This term of the current 
multipole can be a more efficient source of GWs than the mass 
quadrupole, which is relatively small because AP/p is non van- 
ishing only in a narrow region. On the other hand, the integral 
which defines the current multipole extends on a larger vol- 
ume. From simple order of magnitude estimate we found that 
its contribution to the emitted power can be up to one order of 
magnitude larger than the one associated with the mass multi- 
pole. 

4.2.3. Frequency, time duration and amplitude of the 
GW bursts 

Concerning the value of the frequency of the emitted GW 
bursts, in the simplified Newtonian scheme we have adopted, 
the eigenfrequencies of the f-mode and the p-mode are re- 
spectively of 3 and 8 kHz. Calculations of these frequencies 
in full General Relativity, give smaller values, typically of the 
order of 2 kHz for the f-mode ptenhar e t al. (2004) |. The cor- 
responding duration is tj ~ 0.18 s l lBenharTTal. (2004)] ) and 
therefore the total energy emitted by the f-mode is of order 
Ef ~ 1.2 X 1O~ 6 M . 

In the estimate of the f-mode frequency we have up to 
now neglected the effect of rotation. It is known that in a ro- 
tating star the eigenfrequencies are shifted towards lower val- 
ues (Ferrari et al. (2004) I. In particular, the frequency of the f- 
mode reduces to a value of order 1 kHz for fast rotating stars. 

From the emitted power it is possible to estimate 
the initial amplitude ho of the GW using the equation 
( |de Freitas Pacheco (1998p : 

.1/2 



1 



df j Du) f 

where D is the distance of the source. We obtain therefore 



(17) 



h 







3x10-*^ ^ 



lkHzj 



(18) 



The time dependence of the GW is given by 
h{t) = h acp-Qfrr-W , 
and the Fourier transform reads therefore: 

1/t/ 2 + (a> - a>/) 2 

If w - ojf ~ 1000 s- 1 then h ~ (10~ 6 s 2 )/jo/t/. Instead if co - 
(Of ~ 1/t/ then h ~ %r//2. 



(19) 



(20) 



It can be of interest to compare our results with the 
characteristics of the signal detected by NAUTILUS and 
EXPLORER. The experimental value of the Fourier transform 
of the amplitude is h ~ 2 x 10~ 21 s. Clearly our model can ap- 
proach this value only if the frequency <x>f is very near to the 
resonance frequency of the two gravitational bars, which is of 
order of 1 kHz. Even in this case, to obtain in our model an 
amplitude comparable to the experimental one it is necessary 
to assume that the emitted power is significantly larger, maybe 
due to the current multipole term and/or to the effect of using a 
more realistic model. 

Finally, we want to stress that, even using the toy model, 
the energy released in GW bursts is larger than the energy of 
the bursts in other models of GW bursters, as e.g. soft gamma 
repeaters. There, an energy ~ 10~ 10 - 1O~ 9 M can be released, 
but only under the assumption that all the elastic stress is com- 
pletely converted into GW bursts ( |de Freitas Pachec~1 998)l. 



5. Conclusions 

We would like finally to discuss the phenomenological rele- 
vance of the proposed model. Concerning the stellar objects 
which can be possible candidates for our model, if we assume 
that all neutron stars are born with a large value of angular ve- 
locity, they will all enter the instability window as described 
in the first scenario. Taking a neutron star production rate of 
order 0.02 per year in our galaxy, and assuming the possibil- 
ity to detect GW bursts up to a distance of order 1 kpc, the 
probability of finding an active burster in this region is of order 
percent, if the total duration of the emission phase is of order 
50 years. Moreover, if the possibility of "recycling" described 
in the second scenario is taken into account, the probability can 
be larger. A precise estimate of the probability would require a 
precise knowledge of the number of accreting millisecond pul- 
sars, what is not known at the moment. 

An important feature of the model discussed here is that no 
neutrino signal is expected in detectors as LVD (in the Gran 
Sasso National Laboratory) or even in more massive detectors 
as Super Kamiokande. In our model an energy Etrans ~ 1O~ 4 M 
is deposited inside the star, near the metastable layer, during 
the transition. The neutrinos produced by URCA processes 
have typical energies of a few MeV and they will scatter many 
times inside the star, degrading their energy before escaping. 
Therefore the emitted neutrinos have energies below the thresh- 
old of neutrino detectors. 

Finally, let us discuss possible signatures of the mechanism 
here proposed. Concerning the number of epochs during which 
bursts are emitted, we can expect that it will decrease from a 
few active periods in ten years to a few periods in the next hun- 
dred years, because the angular velocity decreases by roughly 
the same amount in the first ten years as in the following hun- 
dred years. The amplitude of the signal should also decrease 
because, as shown in Figs.@]and[5] when the angular velocity 
is lower, the variation of the radius of the star is also smaller. 
Both these features are related to the draining of the bursts en- 
ergy source which, in this model, is ultimately the rotational 
energy of the star. 
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